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Abstract 

A narrow dip structure has been observed at 1 .9 GeV/c^ in 
a study of diffractive photoproduction of the Stt+Stt^ final 
state performed by the Fermilab experiment E687. Prelim- 
inary results on new structures in 27r+27r^ final states are 
also given. 

1 INTRODUCTION 

The Fermilab Experiment 687 collaboration has collected 
a large sample of high-energy photoproduction events, 
recorded with the E687 spectrometer [^[^ during the 
1990/91 fixed-target runs at the Wideband Photon beam- 
line at Fermilab. Although the experiment is focussed on 
charm physics, a very large sample of diffractively photo- 
produced light-meson events was also recorded. This paper 
reports on a study of the diffractive photoproduction of 
the Stt+Stt^ final state and the observation of a narrow dip 
in the mass spectrum at 1.9 GeV/c^. 

Preliminary results on diffractive photoproduction of 
27r"'"27r^ final states are also shown. 

2 E687 SPECTROMETER 

In E687, a forward multiparticle spectrometer is used to 
measure the interactions of high-energy photons on a 4- 
cm-thick Be target. It is a large-aperture, fixed-target spec- 
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trometer with excellent vertexing, particle identification, 
and reconstruction capabilities for photons and tt^'s. The 
photon beam is derived from the Bremsstrahlung of sec- 
ondary electrons of « 300 GeV endpoint energy, which 
were produced by the 800 GeV/c Tevatron proton beam. 
The charged particles emerging from the target are tracked 
by a system of twelve planes of silicon microstrip detec- 
tors arranged in three views. These provide high-resolution 
separation of primary (production) and secondary (charm 
decay or interaction) vertices. The momentum of a charged 
particle is determined by measuring its deflections in two 
analysis magnets of opposite polarity with five stations of 
multiwire proportional chambers. Three multicell thresh- 
old Cerenkov counters are used to discriminate between pi- 
ons, kaons, and protons. Photons and neutral pions are re- 
constructed by electromagnetic (EM) calorimetry. Hadron 
calorimetry and muon detectors provide triggering and ad- 
ditional particle identification. 

3 EVENT SELECTION 
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Figure 1: Distribution of Stt+Stt^ invariant mass after ap- 
plying a cut on the total energy deposited in the calorime- 
ters with respect to the total energy in the spectrometer. 

Pions are produced in photon interactions in the Be tar- 
get. The data acquisition trigger requires a minimum en- 



ergy deposition in the hadron calorimeters located behind 
the electromagnetic calorimeters and at least three charged 
tracks outside the pair region. The microstrip system and 
the forward spectrometer measure the six-pion final state 
(in this paper,67r refers to the Stt+Stt^ state) with a mass 
resolution cr = 10 MeV/c^ at a total invariant mass of 
about 2 GeV/c^. It is required that a single six-prong ver- 
tex be reconstructed in the target region by the microstrip 
detector, with a good confidence level. Such a requirement 
rejects background due to secondary interactions in the tar- 
get. Exclusive final states are selected by also requiring that 
the same number of tracks be reconstructed in the magnetic 
spectrometer The six microstrip tracks and the six spec- 
trometer tracks are required to be linked, with no ambigu- 
ity in associating the microstrip and spectrometer tracks. 
Events with particles identified by the Cerenkov system as 
definite electrons, kaons, or protons, or as kaon/proton am- 
biguous are eliminated and at least four out of six parti- 
cles have to be positively identified as tt*. Particle iden- 
tification is tested by assuming that one or two out of the 
six tracks is a K^, by computing all two-track invariant 
mass combinations, and verifying that there is no evidence 
of a peak at the K* or at the (j> mass. We eliminated fi- 
nal states with tt^'s by rejecting events with visible en- 
ergy in the electromagnetic calorimeters that was not as- 
sociated with the charged tracks. A cut in this variable 
(£'cai/-E'67r < 5%) is applied on the calorimetric neu- 
tral energy normalized to the six-pion energy measured in 
the spectrometer. The distribution of the six-pion invariant 
mass after these cuts is shown in Fig. 1. The plot shows a 
structure at 1.9 GeV/c^. In the following, only the 6tt mass 
region around this structure will be analyzed. 

For diffractive reactions at our energies, the square of 
the four-momentum transfer t can be approximated by the 
square of the total transverse momentum of the diffrac- 
tively produced hadronic final state. Using this definition, 
the distribution of diffractive events. Fig. 2, is well de- 
scribed by two exponentials: a coherent contribution with 
a slope be — 54 ± 2 (GeV/c)~^ consistent with the Be 
form factor and an incoherent contribution with a slope 
b^ = 5.10 ±0.25 (GeV/c)-2. 

Taking only events with P^ < 0.040 GeV^/c^, we eval- 
uated a contamination of about 50% from nondiffractive 
events. This incoherent contribution shows no structure in 
the 1.2—3.0 GeV/c^ mass range. Fig. 3. The diffractive 
mass distribution was obtained by subtracting this contri- 
bution, parametrized by a polynomial fit, and dividing the 
yield by the detection efficiency. 

The detection efficiency was computed by modeling 
diffractive photoproduction of a mass M, using the experi- 
mentally found slope be, and simulating the decay Gtt 
according to phase space There is no threshold or dis- 
continuity for the efficiency, particularly in the region of 
the dip structure. At 1.9 GeV/c^, the (self-normalized rel- 
ative) efficiency A varied as dA/A/dMej^ = 0.15/ GeV/c^. 
The efficiency and the efficiency-corrected distribution of 
the six-pion invariant mass for diffractive events, in the 
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Figure 2: Tranverse momentum squared distribution 
showing the diffractive peak. 



mass range 1.4—2.4 GeV/c , are shown in Fig. 4. There 
no evidence, albeit with large combinatorial backgrounds. 



for resonance substructure, e.g., 



in the 6tt 



data below AfgTi- = 2.0 GeV/c , either at the mass re- 
gion of the dip or in nearby sidebands. Similarly, the ef- 
ficiency or acceptance exhibited no threshold, edge, or dis- 
continuity over the entire mass region observed, when the 
six pion state was simulated as a sequence of decays of in- 
termediate two-body resonances, for example + 
(p°7r+) + (p°7r^) (7r+7r^7r+) + (vr+TT^vr^), even un- 
der extreme assumptions of full longitudinal or transverse 
polarizations for the initial state. 

The presence of a dip at 1.9 GeV/c^ was verified by 
several systematics checks. Diffractive photoproduction of 
pairs or the associated production of plus a 
charm baryon at low t (where the decay products of the 
other charm meson or charm baryon are missed) followed 
by the 6 tt decay of the D" or are estimated by Monte 
Carlo simulation to be negligible contributions. It was also 
checked that demanding more stringent cuts (i.e., requiring 
that all six particles be identified as tt* , applying a sharper 
cut on the calorimeter neutral energy, or subtracting the in- 
coherent contribution bin by bin) increases the statistical 
errors without significantly affecting the behavior shown in 
Fig. I 

A three-parameter polynomial fit was performed, solid 
line in Fig.4, to explore the hypothesis that any structure in 
this distribution is a statistical fluctuation. The normalized 
residual distribution, evaluated for each 10-MeV/c^ bin, is 
good in the full invariant mass range, with the exception 
of the interval centered at 1.9 GeV/c^, the region of the 
claimed dip, where a poor ^ 1Q~^ confidence level inter- 
val, is obtained, making it highly unlikely that the observed 
dip is a statistical fluctuation. Incoherently adding a Breit- 



Wigner to the fit does not improve the fit quality much. 
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Figure 3: Distribution of Stt+Stt^ invariant mass in the 
1.2—3.0 GeV/c^ mass range: coherent plus incoherent 
contribution. Dotted distribution: incoherent contribution. 



FITTING THE SIX-PION INVARIANT 
MASS DISTRIBUTION 
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Figure 4: Acceptance-corrected distribution of Stt+Stt" 
invariant mass for diffractive events after subtracting inco- 
herent contribution. Solid line: second-order polynomial 
fit. Dashed line: polynomial fit with incoherently-added 
Breit-Wigner. Upper dot line: relative detection efficiency 
(arbitrary normaUzation). 

Because of the narrow width, the E687 spectrome- 
ter mass resolution, ct = 10 MeV/c^ at 2 GeV/c^, 



was unfolded by applying the method described in 
Ref. [|[]: the experimentally observed data distribution r(a;) 
and the unfolded mass distribution a{x) are related by 
a{x) ~ r{x) — 0.5(7^ • r{x) , where r{x) is the sec- 
ond derivative with respect to M6,r for the observed distri- 
bution. This relationship results from applying a Fourier 
transform and approximating the resolution function by 
g{x) — exp{— ^1^1 ), which is marginally different from 
a Gaussian. A fit similar to the following one is used to 
obtain r{x) from the unfold ed data. The data after un- 
folding are shown in Fig. ^. 

The dip structure at 1.9 GeV/c^ has been character- 
ized by a two-component fit, adding coherently a rel- 
ativistic Breit-Wigner resonance to a diffractive contin- 
uum contribution. The continuum probabiUty distribution 
Fjs [M] has been modeled after a Jacob-Slansky diffrac- 
tive parameterizationP], plus a constant term cq 



Fjs{M)= fjs{M)^co+ci 



The Jacob-Slansky amplitude fjs{M) is assumed to be 
the purely real i4>js = 0) square root of the probability 
function Fjs{M). For the fit, a relative phase factor e**^, 
independent of mass, and a normalizing factor multi- 
plied a relativistic Breit-Wigner resonance term, giving the 
overall amplitude 



A{M) = fjs{M) + ar 



-MrTe 



AP - M2 + iMrT 



Fit results are shown in Table [l] and in Fig. || for a fitted 
mass range from 1.65 to 2.15 GeV/c^, symmetric with re- 
spect to the dip. Quantities shown are the mass and width 
of the resonance, the amplitude ratio Or/ f^siMi ) between 
the Breit-Wigner function and the Jacob-Slansky contin- 
uum, the relative phase and the x^/dof. Fit values show 
consistent evidence for a narrow resonance at = 1.91 1 
± 0.004 ± 0.001 GeV/c^ with a width F = 29 ± 11 ± 
4 MeV/c^, where the errors quoted are statistical and sys- 
tematic, respectively. The fit values shown in Fig. |] and 
represented by the parameters of Table 1 are stable with 
acceptable x^/dof over varying mass ranges from 1.65 to 
2.3 GeV/c^. We quote as systematic error the sample vari- 
ance of the fit values due to our choice of fit mass range. 
The quality of the fit deteriorates somewhat as the upper 
limit of the fit for this simple model is extended from 2.1 
to 3.0 GeV/c^. However, the only fit parameter that is af- 
fected is the width, which varies from 29± llMeV/c to 
40±20MeV/cl 

5 DIFFRACTIVE PHOTOPRODUCTION 

OF 2(7r+7r-) EVENTS 

About one million of 2(7r+7r~) events have been recorded 
by E687 in the 1990/91 fixed target runs at the Fermi Na- 
tional Accelerator Laboratory.Their analysis is in progress 



Table 1: Fit results for a mass range from 1.65 to 2.15 



Mr (GeV/c^) 


1.911 ±0.004 


r (MeV/c^) 


29 ±11 


ar//js(MO 


0.31 ±0.07 


<t> (deg.) 


62 ±12 


X^/dof 


1.1 




1.5 1,6 1.7 1.8 l.y 2 2.1 2.2 2.3 2.4 



M,,(GeV/c") 

Figure 5: Acceptance-corrected distribution of Stt+Stt^ 
invariant mass for diffractive events. The mass resolution 
has been unfolded. Fit parameters are listed in Table 1 . 



ratio of the residuals to the total yield in 47r is of the same 
order of the p — ui mixing effect in 27r[^, shown in Fig. ^ 
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and in the following only a preliminary presentation is 
given. As in the Stt+Stt" study, exclusive final states 
are selected based on the total number of charged tracks 
seen in the spectrometer and final states with 7r"'s are re- 
jected by requiring no visible energy in the electromag- 
netic calorimeters. The particle identification has been 
done applying the same cuts described in section 3. The 
t-distribution shows a slope of about 60 (GeV/c)~^ , 
in good agrement with the Stt+Stt" sample. Taking only 
events with P| < 0.0625 GeVVc^, our 2(7r+7r") sam- 
ple remains with a contamination of about 25% from non- 
diffractive events. 

The acceptance corrected 2(7r+7r^) mass spectrum is 
shown in Fig. ^top). Although we expect, to first order, that 
the bump in the mass spectrum below 2.0 GeV/c^ is domi- 
nated by the pi lASQ) vector meson, our fit is not good (solid 
line) even when we add 3 P-Wave interfering Breit-Wigner: 
the is 467.0 for 274 degrees of freedom.As the accep- 
tance is almost flat for M^^t^ 1.6 GeV'^/c'^, the deviations 
observed should be of physical origin. The residuals (data 
- fitted values) shown in Fig. ||(bottom) clearly show struc- 
tures with a non negligible statistical significance. Their 
study is in progress. For the time being we can say that the 



Figure 6: Acceptance corrected 2(7r+7r^) mass (top) and 
residual (data - fitted values) distributions (bottom).The 
solid line is the 3 P-Wave Breit-Wigner fit. 



6 CONCLUSIONS 

The diffractive photoproduction of Stt+Stt^ has been stud- 
ied by E687. Evidence has been found for a narrow struc- 
ture near AfgTr = 1.9 GeW/c^. If this dip is characterized 
as the destructive interference of a resonance with the con- 
tinuum background, then the parameters of this resonance 
would be Mr = 1.91 1 ± 0.004 ± 0.001 GeV/c^, with T = 
29 ± 1 1 ± 4 MeV/c^.Such a resonance could be assigned 
the photon quantum numbers (J^'-' = 1 ) and G=H-1, 1=1 
due to the final state multiplicity. 

The possible presence of other structures in the 2(7r+7r~ ) 
final state needs further study. 
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